Linearized Bekenstein Varying Alpha Models 
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We study the simplest class of Bekenstein-type, varying a models, in which the two available 
free functions (potential and gauge kinetic function) are Taylor-expanded up to linear order. Any 
realistic model of this type reduces to a model in this class for a certain time interval around the 
present day. Nevertheless, we show that no such model is consistent with all existing observational 
results. We discuss possible implications of these findings, and in particular clarify the ambiguous 
statement (often found in the literature) that "the Webb results are inconsistent with Oklo" . 



I. INTRODUCTION 

In theories with additional spacetime dimensions 
there are typically many light or massless degrees of free- 
dom, which can give rise to a number of observable cos- 
mological consequences. Noteworthy among these are 
variations of the fundamental couplings y, |^ (with 
the ensuing violations of the Equivalence Principle [EH) 
and contributions to the energy density budget of the 
Universe. 

In recent years there has been a growing body of evi- 
dence for the presence of these two effects. Type la su- 
ernovae j7|, the Cosmic Microwave Background (CMB) 
and lensing data are all consistent with the ex- 
istence of the so-called dark energy component, whose 
gravitational behaviour is very similar to that of a cos- 
mological constant, and which indeed appears to have 
become the dominant component in the energy budget 
of the Universe at a redshift z ^ 1. 

On the other hand there is some (somewhat more con- 
troversial) evidence for the spacetime variation of the 
fine-structure constant a, coming from both quasar ab- 
sorption systems (at redshifts z ~ 1 — 3, fTolllll IT^) and 
the Oklo natural nuclear reactor {z ~ 0.1, [Ijl)- There is 
also a further claim of a varying proton to electron mass 
ratio, also at z ~ 3 0|. While it is conceivable that 
hidden systematic effects are still contaminating some 
of these measurements, an unprecedented effort is being 
made by a number of independent groups and using a 
range of techniques, to search for such variations at var- 
ious key cosmological epochs, which should soon clarify 
the situation. There is also a range of other constraints, 
either local (from atomic clocks) or at low [T^ (from 
Rhenium decay in meteorites) or high redshift |l7lll8lll9| 
(from the CMB and BBN), with much stringent ones 
forthcoming |20ll2l|. 



It goes without saying that from the point of view of 
a fundamental theory there is more than ample freedom 
to allow the dark energy and the varying couplings to 
be due to different degrees of freedom in the theory, and 
even to emerge through different physical mechanisms. 
Nevertheless, it is useful to study the simplest case in 
which the two have a common origin, as this will in prin- 
ciple have the fewest free parameters and can therefore 
be better constrained. 

In what follows we will discuss a very simple toy model 
for this which, although arguably oversimplified from the 
particle physics point of view, has the advantage of hav- 
ing a minimal number of free parameters. The basic idea 
is to consider Bekenstein-type models , and reduce the 
freedom in the two free functions (the potential ¥{(/)) and 
the gauge kinetic function Bp{(l))) by Taylor-expanding 
them around the present day, and retaining only terms 
up to linear order. In fact, we will see that its free pa- 
rameters are in some sense too few, so that the model 
is very tightly constrained by existing observations, and 
indeed ruled out if all of them are correct. Still the model 
can be useful in providing some guidance for the likely 
requirements of successful, fundamental theory inspired 
models. Other interesting analyses of this class of models 
can be found in H H IM El HHH . 

We will start in Sect. II with a brief overview of the 
Bekensein-type models. We then discuss in more detail 
the linearized regime that interests us (Sect. Ill) and dis- 
cuss it in the context of existing observational data (Sect. 
IV). Finally Sect. V summarizes our results and briefly 
discusses further prospects. Throughout this paper we 
shall use fundamental units with fi = c = G = 1. 



II. OVERVIEW OF BEKENSTEIN MODELS 
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The variation of the fine-structure constant in 
Bekenstein-type theories ^ is due to the coupling of 
a scalar field (j) to the electromagnetic field tensor F^^, 
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through a term of the generic form 



, (1) 



where m ~ mp; — 1 and -Bp , known as the gauge kinetic 
function, is the effective dielectric permittivity. This 
should be explicitly specified in a fundamental theory, 
but can be phenomenologically taken to be a free func- 
tion. The fine-structure constant is then given by 



a(0) 



ao 
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and at the present day one has Bp((j)Q/m) — I. 

The inclusion of an interaction term such as Bp{<j))F'^ 
that is non-renormalizable in 4D requires, at the quan- 
tum level, the existence of an ultra-violet momentum cut- 
off. Any particle physics motivated choice of this cutoff 
will destabilize the quintessence potential, since it will 
yield a mass term much larger than the quintessence one 
(recall that typically ~ Hq). This is therefore a fur- 
ther fine-tuning problem, akin to the cosmological con- 
stant one. In our phenomenological approach (and in 
common with all previous work on these models) we will 
ignore this problem, assuming that any mechanism that 
solves the cosmological constant problem will also solve 
this one. (Such a mechanism must exist if the dark en- 
ergy of the universe is indeed quintessence-like.) 

Assuming that the cosmological change in the scalar 
field (f) is small (at least in recent epochs), one can ex- 
pand all couplings around their present-day values, in 
particular 



Bp - 

m 



1 + C-F 



^0,1, 
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corresponding to a variation of a (again relative the 
present day) 



a 



1 
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Given that the classical predictions of the model will be 
independent of the particular choice for the mass m, we 
shall take m = mpi = 1 throughout. In addition to the 
variation of the fine-structure constant, this coupling is 
responsible for an effective non-universality of the gravi- 
tational force, which through Equivalence Principle tests 
leads to the constraint 



ICfI <5x fO-4; 



(5) 



note, for example, that in Bekenstein's original theory 
one has C,p = —2. 

The evolution of the scalar field is then typically of the 
form 



< + 'iH(t> 



dV 



(6) 



where pm is the matter density of the universe. Given 
a complete particle theory, C,p will be specified and it 
should be possible to calculate the coupling of the scalar 
field to matter, Typically the scalar field will only 
evolve significantly during the matter era — a result that 
is well known from the study of scalar-tensor theories 
|29l Isoj l . A non-zero Qp is in principle sufficient to en- 
sure a cosmological variation of 0, driven by the electro- 
magnetic part of the baryon mass density, and hence a 
variation of a. However, the resulting change will typ- 
ically be small if this is the only source. For example, 
in the original Bekenstein model, one can only fit the 
Webb results at the cost of having a huge violation of 
the weak equivalence principle; conversely if one wants to 
satisfy these constraints the typical allowed variation is 
only Aa/a ^ 10^^*^. This constraint can only be evaded 
by saying tha.t (j) couples only (or predominantly) do the 
dark matter |23|. 

In addition to the equivalence principle constraint, 
there are a number of bounds or detections which re- 
strict the cosmological evolution of a, as already dis- 
cussed above. Moreover, if one assumes that 4> is also 
providing the dark energy its evolution will be further 
constrained through its present contribution for the en- 
ergy budget and the evolution of its equation of state 

(siiiiiia. 



III. ANALYSIS OF THE LINEARIZED CASE 

Let us consider the class of models of a neutral scalar 
field coupled to the electromagnetic field with 



where 



C — + L^p + £othcr 



= -9^09^0 - y(0) , 



-<i>F 



(7) 



(8) 
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and Mother is the Lagrangian density of the other fields. 
We will make the simplifying assumption that both V{(j>) 
and a are linear functions of 4>, namely 



vici,) = v{M 



dV 



(0 - M 



and 



da 



ao 



(10) 
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with both dV / d(j) and da/d(j) assumed to be constants. 
One then has 



^ = {B,^ - 1) - 
a ao da 



M- (12) 
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We will further assume that dV/d(j) < 0. Assuming 
that the interpretation of the Webb et al. results as ev- 
idence for a variation of the fine structure constant is 
correct this implies a smaller value of a in the past and 
consequently da/d(p > 0, though if the claimed Oklo de- 
tection is also true there must be oscillations. Let us first 
comment on these assumptions. In any model both V{4>) 
and q;(0) can be taken as linear functions of cf) for some 
limited period of time around today. For how long that 
assumption holds is of course model dependent. In this 
paper we will be considering a particular class of mod- 
els for which this assumption is valid for a considerable 
time, possibly even all the way from the epoch of nucle- 
osynthesis up to the present time. In other words, we are 
effectively testing the validity of this assumption for the 
class of Bekenstein models, assuming the validity of the 
claimed low-redshift detections of a varying fine-structure 
constant. 

In a spatially flat Friedman-Robertson- Walker universe 
the equations of motion are given approximately by 



(13) 
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where 



and 



87r((/.V2 -t- ¥{(/))) 
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where Pco is the present-day critical density. Note that 
the last term can be neglected at the present time but 
it becomes important at early times since the main con- 
tribution to this term comes from baryons whose energy 
density varies as a~^. In summary, in this context, a par- 
ticular model is fully specified by the parameters 0o, Vq, 
dV/d(j) and da/ dcj), in addition to the nominal cosmolog- 
ical parameters. 

Let us start with some order of magnitude constraints 
on the value of 5(j) = </>(2 = 0) — 4'{z = 1), with z ~ 1 
being singled out as the approximate redshift at which 
the matter domination epoch ends (do not confuse this 
particular 5(j) with the generic A0 that has already been 
introduced above) . If our theory is to produce interesting 
variations of the fine structure constant with redshift the 
first term in the right hand side of (|17|l needs to be the 
dominant one. Hence we have 



5V ^ Pco 



(19) 



which implies that Scj) < 10~^. On the other hand if we 
want to have a variation in the fine structure constant of 
Sa/a ^ 10~^ without violating the equivalence principle 
one needs 



Sa S(/) da ^q-g 
a a d(j) 



(20) 



Given that the equivalence principle tests give the con- 
straint 



da 



< 10" 



(21) 



'^V2-y(^) 

(^V2 + y(0) 
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Since any variation of the fine structure constant from 
the epoch of nucleosynthesis onwards is expected to be 
very small [itI IT^ we have also neglected the minor con- 
tribution that such a variation has in the evolution of 
the baryon density (included in i~l,no)- The equation of 
motion for the field is given by 



. + 3if (/) : 



dV 



(17) 



It is crucial to discuss the relative importance of the last 
two terms. We shall assume that the main contribu- 
tion to the last term comes from baryons. Given that 
F^^F^,^ = 2{B'^-E^) < the last two terms in equation 
p7|l have opposite signs. It has been shown that the time 
variations of the fine structure constant induced by the 
last term are two small to ever be observed (if Equiva- 
lence Principle constraints are to be obeyed !2J|). How- 
ever, this does not happen with the first term. Hence, 
in order to have interesting variations of a the first term 
needs to dominate at recent times which will happen if 



(18) 



one has necessarily S 
give a limit on S(j) of 

10" 



> 10 ^. These two constraints 



< 



< 10" 



(22) 



which will be further tightened by future tests of the 
equivalence principle. A related limit on dV/d4> gives: 



2 dV 
Ho<-rr< 



-W-^Hl, 



(23) 



In order for the lower limit on dcf) to be satisfied one needs 
w^{z = 0) — 1 > 10^**. It is also straightforward to verify 
that w — > — 1 very rapidly as we move backwards in time. 

Before we investigate the possible role of the scalar field 
responsible for the variation of the fine structure constant 
as a quintessence candidate it is instructive to study the 
evolution of the scalar field assuming that its contribution 
to the dynamics of the universe is subdominant. From 
the equations above it is easy to show that 



d^<l> d<P 1 



dV 6 



where 



/(a) = ilmOfl ^ + f^Aoa^ + ^rOCL 



(24) 



(25) 
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FIG. 1: The evolution of the value of Aa/a as a function of 
redshift, in re-scaled units with Ca = —tp7^ ~ — Cfttt4^- 

' ^ 



aa 



(26) 



and the last term has been expressed as a function of 
the behaviour of the matter density (with the constant 
Q absorbing the additional parameters). Assuming that 
dV/ d(j) is a constant one can find the following asymp- 
totic solutions in the radiation, matter and A-dominated 
epochs 



— Ar + BrO, 



1 



= AA+BAa- 
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T) f72 
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1 dV 



2 dV 
27 d^' 
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(27) 
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In a 



(28) 



dV 



3D,\oHq [ d(f) ' \y ^ 3 

(29) 

If we neglect the decaying mode it is possible to match 
the solutions deep in the matter era with the solution 
deep in the radiation era era in such a way that both 
and d(j)/da are continuous functions of the redshift. We 
can then find A a, Ba, Am as 




Hq (solid line) 
-5 X 10-" and 



log(1+z) 

FIG. 2: (Top Panel): The evolution of the value of Aa/a 
as a function of redshift for dV/dcf) = — 10"^ 
and dV/d(j> = (dashed line) with ("f ~ 
& — IO^'^Hq. The electromagnetic term gives a negligible 
contribution to the variation of a if \dV/d(j)\ is large enough 
favouring a larger value of a in the past. (Bottom Panel): 
The evolution of the value of Aa/a as a function of redshift 
for dV/d(j> = -IQ-^Jfo- We clearly see that no cosmological 
significant variation of a exists beyond z — 10. 
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dV 
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dV 
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1 + 31n 
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3" 



In- 



mO 



(30) 



3 JIao 
(31) 

As expected |2^ |29|, |30| , we find that the evolution of 
(j) is negligible during the radiation dominated epoch but 
significant during the matter one, and that the onset of 



cosmological constant domination damps this evolution. 
On the other hand, the cosmological evolution in the case 
where the dark energy of the universe is provided by (jj 
itself rather than by a cosmological constant can be ap- 
proximately inferred from the above analysis, since it is 
observationally known that the gravitational behaviour 
of the former must be very close to the latter. 
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FIG. 3; The evolution of the value of + 1 with redshift, z, 
for dV/d4> — —0.35 assuming fJ^o ~ 0.7 and ~ —1. Note 
that + 1 evolves very rapidly toward zero when one moves 
backward in time. 



FIG. 4: Comparison of two typical ex amp les of this class of 
models with observational data: Oklo circles), Rhenium 

decay from meteorites f[T^. vertical bar) and quasar data 
(Murphy et al. [TUl dark shaded box, Chand et al. [l^ light 
shaded box). Either this class of models is not valid up to 
redshifts about z ~ 3, or some of the data is strongly affected 
by systematics. See main text for further discussion. 



IV. DISCUSSION 

In Fig. ^we plot the evolution of the fine structure 
constant as a function of redshift. Given that 4>o is pro- 
portional to dV/d4> and the dynamics of the universe near 
the present time is constrained to be very close to that of 
a universe with fi^, ~ 0.7 and uj^ ^ —1, the shape of the 



evolution curve of Aa/a as a function of redshift, z, is 
unambiguously predicted by our model up to a normal- 
ization factor 



(32) 



We see that for a value of ^ —10^^ (take for example 
— —10"'* and dV/d<j) — — 10"^i?Q one has interesting 
variations of a which we would expect to be able to detect 
some time in the near future. However, note that the 
evolution of a with redshift is still quite significant at 
very low redshifts, which indicates that it is not possible 
to reconcile the Oklo ^3 oi' nieteorite 0| results (at z = 
0.14 and z = 0.45 respectively) with the Webb/Murphy 
results 0, ^3 in the context of our model. If we take 
the Oklo and meteorite limits seriously the maximum 
variation of a that is allowed in this class of models by 
a redshift z ^ 2 say is about 10~^, clearly below the 
Webb/Murphy results, though perhaps compatible with 
Chand et al. [ij. More on this below. 

Examples of two particular models are displayed in Fig. 
121 In the top panel we plot the evolution of the value of 
Aa/a as a function of redshift for dV/dcj) — —IO^^Hq 
(solid line) and dV/d(() = (dashed line) with C,f — 
—5 X 10~* and 6 — IO^^Hq. For a nearly flat scalar 
field potential the electromagnetic term in equation H17|) 
is the main source of a variation of the fine structure con- 
stant, favouring a larger value of a in the past. However, 
in that case the variations are too small to be of any cos- 
mological significance. In the bottom panel we plot the 
evolution of the value of Aa/a as a function of redshift 
for dV/d(j) = —10~^Hq showing that no cosmological sig- 
nificant variation of a beyond z — 10. 

Given a value of Hq, the parameters Vq and 0o deter- 
mine in a unique way w^i^q and fi^o- Hence, the equation 
of state Wci,{z) is a function of dV/d(f>, fJ^o ^^id lu^q only 
and will evolve very rapidly towards a cosmological con- 
stant with Lj^ —^ —1 when we move backwards in time. 
This is clearly shown in Fig. |31 (again, note that w^o 
and ri^o E^re constrained to be very close to —1 and 0.7 
respectively) . 

Let us now go back to the issue of the comparison with 
observational data, which is plotted in Fig. 0] against 
two typical models. Note that the Oklo |13] and me- 
teorite 01 data are plotted as points, since they apply 
to specific redshifts. However, for the quasar absorp- 
tion data 0, we've chosen to plot them as bands, 
rather than plotting the individual observational points 
(or some binning thereof). This choice is partially moti- 
vated by the fact that the error bars are still much larger 
than those from Oklo and the meteorites, and also be- 
cause the quasar 'distilled' results are often quoted as a 
single number that is supposed to apply to a range of 
redshifts. Note however that this practice can be mis- 
leading. For example, a number of authors quote the 
latest results by Webb and collaborators jTJ as 



Aa 



a 



(-0.54 ±0.12) X 10" 



0.2 < z < 3.7, 



(33) 
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ignoring the fact that even though the sample spans all 
that redshift range, the data only prefers a value of a 
different from today's beyond redshift z ~ 1 — hence our 
choice for the horizontal range of the (dark grey) band. 
For the Chand et al. results (light grey bands) we have 
included their two possible results, 

^ = (-0.06 ±0.06) X 10"^ 0.4<z<2.3, (34) 
a 

assuming terrestrial isotopic abundances (case 1), or 

— = (-0.36 ±0.06) X 10"^ 0.4<z<2.3, (35) 
a 

assuming low-metalicity isotopic abundances (case 2). 
The true result is expected to lie somewhere between the 
two. 

A number of interesting points can be inferred from 
the figure. As a first remark, let us point out that even 
though the Webb/Murphy and the Chand et al. results 
seem to be statistically inconsistent, the discrepancy may 
be much smaller than one could guess by simply compar- 
ing (|33|) and H34|l. Nevertheless, the question remains as 
to whether there are hidden systematics contaminating 
one or both of the data or the analysis pipelines. 

We have already pointed out that the distinguishing 
feature of this class of models is that significant varia- 
tions of a occur relatively near the present epoch. (As 
a side remark we note that in this context this justifies 
the commonly used assumption of a uniform value of a 
throughout the last scattering epoch when constraining 
variation of the fine structure constant with CMB obser- 
vations 

EmmilHiJl-) However, this late variation 
has dramatic consequences. Roughly speaking, depend- 
ing on the model parameters one can divide models in 
this class into two different types: they can either be 
consistent with Oklo+meteorites+ ()34ll but be inconsis- 
tent with (|33|l . or else be consistent with 133|l -t- H35|l . but 
inconsistent with Oklo+meteorites. 

In other words, if we assume that our linearized class 
of models holds true up to at least redshift z ~ 3 or so, 
then the Webb/Murphy results are indeed inconsistent 
with Oklo/meteorites. We note that a number of authors 
have in the past made the unqualified statement that "the 
Webb results are inconsistent with Oklo" . Such a state- 
ment is not meaningful per se, since any such comparison 
is necessarily mo del- dependent: one needs to specify a 
timescale as well as a model for the redshift evolution of 
a (see the discussion in 0). Indeed one can build models 
where the two can be made compatible — an example is 
[26|. Having said that, here we do find that the two are 
inconsistent for the models we considered. 

This therefore calls for improvements on the existing 
observational results. Following the controversy gener- 
ated by the quasar data results, at least five (to our 
knowledge) independent groups are currently working on 
the subject, using a variety of different methods, so there 
is hope that the situation will be clarified soon. No simi- 
lar interest exists for the Oklo or meteorite data, though 



independent confirmation of both of these results would 
be much welcome since as we have seen they are quite 
more constraining, particularly for the class of models 
that we have discussed. We note that measurements of 
a using quasar data are, notwithstanding the possible 
sources of observational systematics, quite straightfor- 
ward in the sense that one measures a directly. On the 
other hand, measurements using Oklo and meteorite data 
are indirect: what one measures directly here is some 
combination of various couplings, and using them to ob- 
tain constraints on a requires either assuming that other 
couplings don't vary (which is almost certainly unrealis- 
tic) or assuming some (necessarily model-dependent) re- 
lations between them. It is therefore important to check 
how robust these constraints are to the specific assump- 
tions being made to obtain them 

Of course, if both of these observational results survive 
further scrutiny, then our toy model cannot be correct. 
We emphasize again that any realistic model will reduce 
to a model in this class for some period of time close to 
today, so that would indicate that our linearized approx- 
imation will break down very close to today, arguably 
much earlier than one would have thought. Note also 
that this class of models, with a linearized behaviour for 
the scalar field, are arguably the simplest possible mod- 
els for a varying a. Certainly models where a has a 
linear dependence on redshift or on cosmic time (which 
have been explicitly or implicitly assumed by a number 
of authors) are much more un- natural, and it's hard to 
see how such could be obtained from a sensible particle 
physics theory an a way that would be consistent with 
other observational and experimental constraints. 



V. CONCLUSIONS 

We have studied the simplest class of Bekenstein-type, 
varying a models, and compared them to existing obser- 
vational constraints. These are models in which the two 
available free functions (the potential and the gauge ki- 
netic function) are Taylor-expanded around present-day 
values, with terms kept only up to linear order. Despite 
their apparent simplicity, they are interesting to the ex- 
tent that any realistic model of this type should reduce to 
a model in this class for a certain time interval around the 
present day. Nevertheless, their simplicity means that 
very specific predictions ensue, that can be compared 
with existing data. We have shown that no such model 
is consistent with all the existing observational results. 
Hence either some of these observations are dominated by 
unknown systematics or our linearity assumption breaks 
down on a timescale significantly smaller than a Hubble 
time. 

Given that a scalar field that produces a varying fine- 
structure constant can also make a significant contribu- 
tion towards the dark energy of the universe, it's inter- 
esting to speculate on the possible relation between the 
above observation and hints for a time-varying equation 



7 



of state of dark energy.Indeed in the latter context it has 
been argued that something analogous seems to happen: 
observational data seem to disfavour not only a constant 
equation of state, but even a mildly varying one, say 
with a linear dependence in redshift [s^ ,3^ ,3^J. It is 



unclear if the two things are somehow related, but it has 
been said that a coincidence is always worth noticing — 
one can always discard it later if it turns out to be just 
a coincidence. 
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